####################################
# Main effect RDlocrand
####################################

rm(list=ls())

library(Hmisc)
library(ggplot2)
library(stargazer)
library(foreign)
library(rdrobust)
library(rdd)
library(rdlocrand)
library(readstata13)

################
# Prepare data 
################

# read data
load("~/Dropbox/Gender Chile/08_replication/gender_chile_2020aug28.RData")   
names(d)

# Select predetermined covariates to be used for window selector
X = cbind(d$pinera06_2_p,
          d$bachelet06_2_p,
          d$idh98,
          d$rankingsalud98,
          d$rankingeducacion98,
          d$rankingingreso98,
          d$countysize,
          d$distance)

# Running variable and outcome variable
R = d$margin
Y1 = d$permanent_temporary
Y2 = d$female_permanent_temporary
Y3 = d$male_permanent_temporary
Y4 = d$sharefemale_permanent_temporary
D = as.numeric(R>=0)

#########################################
## Select bandwidths
#########################################

set.seed(123456789)

## Window selection with default options
tmp = rdwinselect(R,X)

tmp <- rdwinselect(R,X,wmin=.01,wstep=.0025,reps=10000)

#####################################################################
## Randomization inference Window 1
#####################################################################

tmp = rdrandinf(Y1,R,wl=-.010,wr=.010)
tmp = rdrandinf(Y2,R,wl=-.010,wr=.010)
tmp = rdrandinf(Y3,R,wl=-.010,wr=.010)
tmp = rdrandinf(Y4,R,wl=-.010,wr=.010)

#####################################################################
## Randomization inference Window 2
#####################################################################

tmp = rdrandinf(Y1,R,wl=-.012,wr=.012)
tmp = rdrandinf(Y2,R,wl=-.012,wr=.012)
tmp = rdrandinf(Y3,R,wl=-.012,wr=.012)
tmp = rdrandinf(Y4,R,wl=-.012,wr=.012)

#####################################################################
## Randomization inference Window 3
#####################################################################

tmp = rdrandinf(Y1,R,wl=-.015,wr=.015)
tmp = rdrandinf(Y2,R,wl=-.015,wr=.015)
tmp = rdrandinf(Y3,R,wl=-.015,wr=.015)
tmp = rdrandinf(Y4,R,wl=-.015,wr=.015)

#####################################################################
## Randomization inference Window 4
#####################################################################

tmp = rdrandinf(Y1,R,wl=-.018,wr=.018)
tmp = rdrandinf(Y2,R,wl=-.018,wr=.018)
tmp = rdrandinf(Y3,R,wl=-.018,wr=.018)
tmp = rdrandinf(Y4,R,wl=-.018,wr=.018)
